Floating and sinking: the imprint of massive scalars around rotating black holes 
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We study the coupling of massive scalar fields to matter in orbit around rotating black holes. 
It is generally expected that orbiting bodies will lose energy in gravitational waves, slowly inspi- 
ralling into the black hole. Instead, we show that the coupling of the field to matter leads to a 
surprising effect: because of superradiance, matter can hover into "floating orbits" for which the 
net gravitational energy loss at infinity is entirely provided by the black hole's rotational energy. 
Orbiting bodies remain floating until they extract sufficient angular momentum from the black hole, 
or until perturbations or nonlinear effects disrupt the orbit. For slowly rotating and nonrotating 
black holes floating orbits are unlikely to exist, but resonances at orbital frequencies corresponding 
to quasibound states of the scalar field can speed up the inspiral, so that the orbiting body "sinks" . 
These effects could be a smoking gun of deviations from general relativity. 



I. Introduction. Massive scalars are ubiquitous in 
physics. For example, light scalars spanning several or- 
ders of magnitude in mass are predicted in string-theory 
scenarios Massive scalars are observationally vi- 

able in scalar-tensor generalizations of Einstein's general 
relativity [|| and can be regarded as an effective propa- 
gating degree of freedom in f(R) theories @, @- In this 
paper we consider generic massive scalar fields coupled 
to matter in orbit around a rotating black hole (BH). 

A well-known phenomenon in BH physics is the Pen- 
rose process (for particles) and the associated superradi- 
ant amplification (for waves) @, H| ■ Consider a Kerr BH 
of mass M , angular momentum J = aM and horizon ra- 
dius r + = M + \/M 2 — a 2 , so that the angular velocity 
of the horizon = a/2Mr + (here and below we set 
G = c = 1). A wave with frequency u < uiQh incident 
on the BH (where m is the azimuthal quantum number) is 
amplified in a scattering process, the excess energy com- 
ing from the BH's rotational energy. Superradiance is 




responsible for many interesting effects [9l4l6|. Here we 
explore the interesting possibility that an object in or- 
bit around a rotating BH may excite superradiant modes 
to appreciable amplitudes. As the object orbits around 
the BH it loses energy in gravitational waves, slowly spi- 
ralling in, as shown experimentally by the Hulse- Taylor 
binary pulsar. This follows from energy balance: if the 
orbital energy of the particle is E p , and the total (gravi- 
tational plus scalar) energy flux is Et — E 9 + E s , then 



E v + E 9 + E s 



0. 



(1) 



Usually E 9 + E s > 0, and therefore the orbit shrinks 
with time. However it is possible that, due to superradi- 
ance, E 9 + E s = 0. In this case E v = 0, and the particle 
can hover in a "floating orbit" [1, [l(| . Here we show that 
floating orbits, for which the net gravitational energy loss 
at infinity is entirely provided by the BH's rotational en- 
ergy, can exist for a wide range of scalar-field masses. 
Orbiting bodies will float until they extract sufficient an- 
gular momentum from the BH or until disruptive (per- 
haps nonlinear) effects stop the process. When the BH 
rotates slowly the condition for superradiance at these 
resonances is not met, but we show that resonances at 
small orbital frequencies (corresponding to large positive 
scalar fluxes going into the horizon) still exist, and that 
they cause the object to inspiral faster. 
II A. Setup. The process we consider is quite general. 
It occurs in all theories of gravity with Kerr BHs as back- 
ground solutions and a scalar field of mass hfi s coupled to 
matter (see e.g. [13, EH)- At first order in perturbation 
theory, the field equations for the scalar field reduce to 



FIG. I. Pictorial description of floating orbits. An orbit- 
ing body excites superradiant scalar modes close to the BH 
horizon. Since the scalar field is massive, the flux at infinity 
consists solely of gravitational radiation. 
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aT. 



(2) 



Our main results will be to a large extent independent of 
the source term on the right-hand side, but for concrete- 
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ness we focus on source terms of the form 
df 



T = 



p 5^ (x - X{f)) , (3) 



corresponding to the trace of the stress-energy tensor of a 
point particle with mass m p , where is the background 
(Kerr) metric. In scalar-tensor theories, for example, 
a = 1/871-/ (2 + wbd) (s — 1/2), where wbd is the Brans- 
Dicke (BD) parameter^ and s is an object-dependent 
"sensitivity factor" 0, El ■ 

Weak-field gravitational radiation circularizes the or- 
bit (see below for a proof in the present context), so we 
consider equatorial circular orbits around a Kerr BH, 
but most results apply to more general orbits. Using 
the "adiabatic approximation" we assume that the ra- 
diation reaction timescale is much longer than the or- 
bital timcscalc, and compute the total energy flux Et 
for geodesic orbits. For prograde orbits, energy, angular 
momentum and frequency of a particle at r = ro read 



E„ 



aVM + V^(r - 2M) 



3/ 4 



m p , 



2aVM + ^(r - 3M) 



M (r 2 - 2a^/TTr~^+a 2 ) 
r 3 /4 y/2aVM + y/^(ro-3M) 



m p , 



M 



(4) 
(5) 
(6) 



The four- velocity of the particle on a timelikc geodesic 
reads r 2 m p U a = {{rl + a 2 )Q / A + a(L p - aE p ), 0,0, L p - 



Q = ( 



aE p + aQ/A), where A = r 2 - 2Mr 
a 2 )E p — aL p . 

IIB. Wave emission. Because of the coupling to 
matter, the orbiting object emits both gravitational and 
scalar radiation. Gravitational radiation can be com- 
puted using Teukolsky's formalism [22|. The relevant 
equations and their solution are presented by Detweiler 
[23[ . Here we focus on scalar wave emission. Defining 

<p(t, r, O p ) = W due^-^ Xlm ^ r) 2 S lm {0) , (7) 



we get the non-homogeneous equation for the scalar field 



dr 2 



V 



1 Measurements of the Shapiro time delay require ojbd > 40, 000 
for fi B = [2(1 . but couplings of order lubd ~ O(l) are obser- 
vationally allowed when /i s >10 — 17 eV, and no bounds on ojbo 
exist when [i s ^ 10~ 16 eV 12111 . Considering a supermassive BH 
of mass M ~ 1O 5 M0 and a typical sensitivity s ~ 0.2, these 
bounds translate into a < 8 ■ 10 -3 when fj, s = 0, a < 0.9 when 
fj, s M > 10 -2 and no bounds on a when fi s M > 0.1. 



where dr/dr* = A/(r 2 + a 2 ), 

(X 

Timuj = -j^S lm (n/2)5(r - r ) m p S{mQ p - u) , (9) 

and the effective potential for wave propagation V is 
given (e.g.) in [19j . Let us consider two independent 
solutions X, * and X^ u to the homogeneous equation 
satisfying the following boundary conditions: 



Imco 



as 



00, r+ , 



where kji — to — itiQh and koo = y^oj 2 — fi 2 . Let W 
be their Wronskian. The fluxes of scalar energy at the 
horizon and at infinity are 



E'=mn p k H , 00 \^\ 2 .. 



(10) 



Irnuj 



-a 



X ^(ro)S L {,/2) mp/M ^ (u) 



WW ^/rj+a 2 



IIC. Analytic solution at low frequencies. The 

scalar flux at infinity can be computed in the low- 
frequency regime [2~H ] . For rp /M 3> 1 and I = m = 1, 



Et 



, 2 M 2 (1 - uy /M) 
12tt r 4 



3/2 



-m 2 9(r!p-M s ), (12) 



where O(x) is the Heaviside function. For generic modes, 
at large distances and for uj = mfl p > fi s , scalar ra- 
diation dominates over gravitational radiation: compare 
Eq. (fT2"|) with the standard quadrupole formula E^ = 
32/5 {r /M)~ 5 ml/M 2 . This result is oblivious to the 
presence of the rotating BH. In fact, for u) > u s , the 
fluxes at the horizon are negligible. However, for fre- 
quencies close to n s , a resonance occurs at (25j : 



Ms - Ms 



fx s M 
l + l + n 



, n = 0, 1,... 



(13) 



From Eq. (fT0|) we see that £y < in the superradiant 
regime (fc# < 0). Close to resonance we get (cf. also (2f|) 



A' 



Imuj 



r l+l e -nMr/(l + l+n) 



(14) 



We have verified this result numerically, finding very 
good agreement with the analytical prediction. For the 
fundamental mode n = 0, at resonance, we find: 



W 



— - l+iP T[l + l]T[l + l-2iP] 

+ + l + ' T[2l + l]r[l - 2iP] 



where P = —2Mr + kn / {r+ — r—). Finally we can esti- 
mate the peak flux close to the resonant frequencies. At 
large distances and for I = m = 1, n = we find 



E s r+ 



3a 2 



^m 2 M 



lQ^ + {M 2 -a 2 )[^ + -{^y/ 2 )F 



(15) 



with F = 1 + AP 2 . Quite surprisingly the scalar flux at 
the horizon grows in magnitude with tq and it is negative, 
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ro/M (resonance) 


(am p /M)- 2 E*f™ k 


Q 


: lit 


10- 1 


4.33400288873563 


-0.1828 


1.1 ■ 




10~ 2 


21.4020987080510 


-0.4881 


1.6 ■ 


10~ 3 


10 -3 


99.9339974413005 


-1.1588 


2.3 • 


10" 5 



TABLE I. Orbital radius at resonance and peak scalar flux for 
n = 0, I = m = 1, a = 0.99A4" and several values of fj, s M <C 1. 
For comparison, a typical extreme mass ratio inspiral becomes 
detectable by space-based interferometers at radii ro/M ~ 
50[(10 6 M o /M)(/ cut /10~ 4 Hz)]~ 2/3 , where / cut is the lower 
cutoff for the sensitivity threshold of the interferometer. A 
floating orbit occurs for a > a C rit- Notice that a cr j t is well 
below current observational bounds [2lJ for any /i s . 



due to superradiance, at sufficiently large distances (for 
generic I, the peak flux would scale as Ep peak cx r 21 3//2 ). 
For very small a the peak flux at resonance is instead pos- 
itive, and it can also be very large: for the Schwarzschild 
geometry, 327rA/ 4 ^;P oak - 3a 2 r 2 m 2 . 

III. Floating orbits. From the previous discussion it 
follows that, for any /x s M -C 1, there exists a frequency 
^rcs < Ms fo r which the total flux E^ + £7* + E^ + 
E-ji = 0, because the negative scalar flux at the horizon 
is (in modulus) large enough to compensate for the other 
positive contributions. This expectation is confirmed by 
a full numerical integration of Tcukolsky's equation: sec 
Fig. [2] and Tabic fl] The width of the peak is proportional 
to the imaginary part of the resonant mode u>i oc /x 4i+5 
[25l ] . For 1 = 1, more explicitly, 



24 



(a/M -2n s r+) 



As fi s — > the imaginary part becomes tiny, and an accu- 
rate fine-tuning is needed to numerically resolve the res- 
onance. For example, to resolve the peak at r$ ~ 100A7, 
corresponding to /i s M = 10~ 3 , we tuned the location of 
ro to 25 decimal places. Computing the imaginary part of 
the unstable modes when ^jl s — > is also challenging, but 
we were able to obtain stable results for the resonance 
location and for the height of the peak. A fit to numer- 
ical results for 10M < r < 100M (cf. Table © yields 

V2 • 



tps, peak ^ 



,0.51 

' 



to be compared with ij^P cak ^ 



Eq. (fT5)) . Close to floating orbits 



dE, 



(F F ) dEr 
— ~-(E p -E f ) — 



(17) 



E B =E, 



where Ef is the energy of the particle at the floating 
orbit, and we used the balance condition ([T]). During in- 
spiral, right before reaching the floating orbit, the time 
needed for the particle to increase its binding energy from 
\Ef\ — e to \Ef\ diverges logarithmically. Therefore, float- 
ing orbits are expected to last much longer than a typical 
inspiral timescale, with a potentially striking observa- 
tional signature in the gravitational-wave spectrum. 




FIG. 2. Dominant fluxes of scalar and gravitational energy 
(I — m — 1 and I = m = 2, respectively) for [i a M = 10 -2 , a = 
10 -2 and a = 0.99Af. The inset is a zoom around resonance. 



(16) Sr 



In the adiabatic approximation, the mass and angular 
momentum of the background spacetime are constant. 
However, the negative energy flux at the horizon re- 
duces the BH mass and angular momentum (6M < 0, 
8 J < 0). In order to estimate how long a particle can 
stay in a floating orbit we must go beyond the adiabatic 
approximation. Under ideal conditions, floating would 
stop only when the peak of the scalar flux at the hori- 
zon is too small to compensate for the gravitational flux, 
\E 9 \ > |£p Cak |. From the balance condition ([T]) we find 
that 5E p = 0, which, using Eq. Q, can be written as 

SMr [a 2 + 3(2M - r Q )r + 2ay^/M(r Q - 3M) 

M (3a 2 - 8a/Mro" + (6M - r )r ) 

where we used the relation SM = il p SJ, valid for 
circular orbits. Substituting the equation above into 
Eq. (fl"5j) and approximating SM/8t ~ E* = — E^ = 

-32/5 (r /M)~ 5 m 2 p /M 2 at resonance, we obtain, for 
I = m = 1 and in the limits fx s — > and a 3> 1 / /i. s , 



ARS.peak 
It 



12a 2 M /m„\ 4 ( M 



5tt a 2 \M 



ro 



which is negative: BH mass loss decreases the height of 
the peak on a timescale 



rps, peak 



5E; 



, peak 



5a / M 
32 lm~ 



(-) 



7/2 



j r+ /St v \ ""p, 

which does not depend on the coupling constant a. 

The delayed inspiral may have observational conse- 
quences. In particular, notice that in the absence of 
scalar fields the evolution of orbital frequency scales as 

tl GR /n 2 p ~ 96/5(m p /M)(A/rg 5 / 3 . (18) 

Close to a floating orbit we find instead that 

fWi„g/tt 2 ~ 32(m p /M) 2 (A/r! p ) 7 / 3 . (19) 
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IIIB. Perturbations of floating orbits. It might be 
anticipated that floating orbits could be very sensitive to 
small perturbations, because the resonance has a small 
width proportional to ujj, which is given by Eq. (|16p . 
One important class of perturbations consists of adding a 
small eccentricity e to the motion. As pointed out in |26| , 
eccentricity will produce frequency sidebands. Because 
these sidebands arc typically far away from resonance, 
scalar emission in these channels can be neglected, and 
we find that e tends to decrease close to floating orbits: 



e 352m p /MW 991 M 2 
~e~ lW UJ V 308^ + (o ' 



(20) 



i.e. circular floating orbits stay circular. The argument 
breaks down at small separations, where a better under- 
standing of eccentricity evolution would be desirable. 

Perhaps the most dangerous perturbation would come 
from small "kicks" taking the object from a trajectory 
with frequency VL p to Q p +SQ p . If 5tt p 3> uji (this is a con- 
servative estimate, as floating itself extends over a larger 
portion of the frequency window) this might stop float- 
ing and superradiant energy extraction. When translated 
into perturbations of the position of the orbiting body, 
which we take to be of size R p , floating would stop for 

Sr a /Rp > 10~ 5 10 dfl p (tom) • For typical extreme- 
mass ratio inspirals of interest for space-based detectors 
M/R p ~ 10 6 and r ~ IQM, so Sr /R p > 10" 5 . 

This rough estimate suggests that when ro/M 1, 
floating is easily stopped by small external perturbations. 
At large distance, even thermal fluctuations of the orbit- 
ing body may disrupt a floating orbit. Note however that 
this estimate is overly pessimistic, as it considers only the 
fundamental mode. For the first few overtone numbers 
n the scalar flux at peak is roughly independent of n. If 
this behavior holds at large n, then the effective allowed 
window SD, p would be larger by powers of l/u 2 s . Further 
investigation of this issue is necessary to understand the 
relevance of floating orbits. 

The small orbiting body was assumed throughout to 
be point-like. In reality there will be tidal interactions 
with the field of the massive BH, which will alter the 
gravitational waveform. For extreme mass-ratios such as 
the one we consider, this effect is negligible [27j. Also, 
tidal effects do not in any way affect the existence of 
floating orbits, but they may have an impact in the dis- 
ruption of floating. It is also conceivable that, on such 
large timescales, nonlinear mode coupling could channel 
energy to higher-order modes, slowly washing out the 
peak in the scalar flux. An analysis of tidal effects and of 
the nonlinear regime is beyond the scope of this paper. 
IV. Sinking orbits. For small rotation, long-lived 
modes of frequency w ~ /i s can be excited outside the in- 
nermost stable circular orbit (ISCO) if fin < oj < ^isco- 
This can happen for < a < 0.36M, and the corre- 
sponding resonances are now stable (cf. the sign change 
in Eqs. (TT5)) and (JTHJ) for u s > Cljj)- For these reso- 
nances, the scalar flux at the horizon is large and posi- 
tive. For tq 3> M (or when ji s AI <C 1) the scalar flux 



dominates over the gravitational flux. From Eq. (fT5|) , us 
ing 5E p /5t = —Et ~ — Ep pcak and neglecting mass anc 
angular momentum loss, we get (for I = m = 1, n = 0) 



sinking 



9a 2 m p /M 



8tt(1 - a 2 /M 2 )(2r + n p - a/M)F 



These resonances are not superradiant in origin. In fact, 
they are present even for nonrotating BHs. For very small 
a, Op inV ing/Op ~ 18a 2 TO p (Mf2j,)~ 3 /(256Af) close to res- 
onance: compare the inspiral rate in general relativity 
without scalar fields, Eq. (TTS)) . Thus, close to stable res- 
onances the orbiting body inspirals much faster, and the 
orbit sinks. Although this effect seems huge at small 
frequencies, its timescale is extremely small, since the 
resonance width ~ uji (cf. Eq. (fTo) 1 ). Indeed, we get 



ir (M 2 - a 2 )F 
27 m„ 



2r+ ^ Ti 



a* 



Note that t s = for a normal mode, at tt p ~ fj, a = &h- 
V. Conclusions. We described an extreme form of en- 
ergy extraction from a Kerr BH when a massive scalar 
field is coupled to a point particle in circular orbit around 
the BH. This is the first example of a phenomenon pro- 
duced by a resonance between orbital frequencies and 
proper oscillation frequencies of the BH. Our results 
apply in principle to stationary background geometries 
different from the Kerr solution. The existence of an 
crgorcgion is mandatory. There are few experimental 
bounds on fx s . Because of supcrradiancc, Kerr BHs 
are unstable under massive scalar perturbations, so (in 
principle) the observation of fast-spinning astrophysical 
BHs can exclude the existence of very light scalars (ax- 
ions) dj]. The strongest instability was found to occur 
when n s = 5.63 x 10" 11 (M/M ) _1 eV on a timescale 
t = 32.8(M/M )s [H. Our results apply to generic 
theories of gravity, a particular instance being BD the- 
ory. Floating may occur and be observed at radii ro < 
50M (cf. caption of Table I). From Eq. (fT5)) . imposing 

£; S ,poak > fig &t ^ = §q M ^ WQ find thfi sma U cst vamc Q f 

a which allows for observable effects. This translates into 
wbd < 10 8 (2 • 10 8 ) for a/M = 0.99 (a/A/ = 0.5). Thus, 
floating orbits could in principle provide much more strin- 
gent constraints on BD theory than those coming from 
Solar System observations. 

Current searches for gravitational waves from coalesc- 
ing compact binaries use post-Newtonian waveform tem- 
plates, and are strongly biased towards general relativ- 
ity. If light scalar degrees of freedom couple to matter, 
binaries may merge in a much more interesting way, and 
current searches based on matched-filtering techniques 
may underperform. In particular, gravitational wave- 
forms would carry a clear signature of floating orbits: 
compare Eq. (TT9")) with the standard general relativistic 
prediction of Eq. (fT8| . Future space-based and advanced 
Earth-based detectors may have the potential to reveal 
the existence of floating orbits and of other surprising 
"anomalies" in the orbital evolution of compact binaries. 
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